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We prove the Lee-Huang- Yang formula for the ground state energy of the 3D Bose gas with repul- 
sive interactions described by the exponential function, in a simultaneous limit of weak coupling and 
high density. In particular, we show that the Bogoliubov approximation is exact in an appropriate 
parameter regime, as far as the ground state energy is concerned. 



I. INTRODUCTION 

We consider a three dimensional system of N interacting bosons in a cubic (periodic) box A 
of side length L, described by the Hamiltonian: 

N 

Hn = -J2^ + ^ £ fflote-x,). (1.1) 

z=l l<i<j<N 

Here G A, i — 1, . . . , N, are the positions of the particles, and Aj denotes the Laplacian with 
respect to Xj. Units are chosen such that fi = 2m = 1, where m is the mass of the particles. The 
interaction potential is taken to be vr (x) = X^nez 3 e~' x+nL '/- R °, and a and R are positive 
constants. The Hamiltonian (1.1) operates on symmetric wave functions in the Hilbert space 
L 2 (A N , dxi ■ ■ ■ dxjy), as is appropriate for bosons. 

We are interested in the ground state energy Eq(N) of (1.1) in the thermodynamic limit when 
N and |A| tend to infinity with the density p = N/\K\ fixed, and in a weak coupling regime 
ao <C min{p -1 / 3 , Rq}. The constant a>o is the first Born approximation to the scattering length 
a of the potential (ao/i^e - !*!/^ , which is defined as usual as a — limixi^oo |x|(l — ^o(x)), 
with ipo a solution to the zero energy scattering equation 

-2A^(x) + He-M/^x) = (1.2) 

Rq 

with boundary condition limix!—^ ip(x) = 1. It is well known that, if ao/i?o <C 1, a/ao can be 
written in terms of a convergent series in powers of a,Q /Rq (Born series) , which will be denoted 
by a = clq + X)fc>i a ki an d whose first non-trivial term is given by 

1 f Kk) 2 5^ al 
ai = -l^/ B 3 dk n^-16iV (L3) 

where 

Kk) = £0 r dx e -| X |/« 0e - lk x = — 8^0 (1 4) 

RoU [l + (ki? )2] 2 

The current understanding of the properties of the ground state of (1.1) is based on the 
pioneering work of Bogoliubov [1], who developed an approximate theory of the ground state 
of weakly repulsive bosons. In the regime 1 3> a/R 3> \J pa 3 3> (a/i? ) 2 , Bogoliubov's theory 
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predicts [7] that the ground state energy Eq(N) of (1.1) in the thermodynamic limit N, |A| — > oo, 
with p = N/\A\ fixed, satisfies 

Km ^D.=Anpa(l + -^Vp^ + o(^)) . (1.5) 

This formula was first derived in [2] and it is known as the Lee-Huang- Yang formula. Our 
goal is to prove that the expression (1.5) is asymptotically correct in a regime such that 
a <C p -1 / 3 <C Ro, that is a weak coupling and high density regime. We shall prove the following 
theorem. 

Theorem 1. Let Y = pa r . There exists a positive constant do, which can be chosen to be 
d a = 1/69, such that, if < d < d a and a/Ro = 0(Y 1 / 2 ~ d ), then (1.5) is valid, asymptotically 
as Y -> 0. 

This result represents the first rigorous proof of the Lee-Huang- Yang formula for the ground 
state energy of a weak-coupling Bose gas. We note that for d < 1/6, R^p 3> 1 and hence 
Theorem 1 concerns the high density regime. Our result is not expected to be optimal. In fact, 
the formula (1.5) is expected to hold even for d = 1/2, i.e., for a/Ro fixed and pa 3 — > [2], but 
the Bogoliubov approximation is not valid in this case. The prediction of Bogoliubov's theory 
is that (1.5) should be valid for any < d < 1/4, i.e., in the regime a/Ro 3> \J pa? » (a/Ro) 2 . 
The latter condition is necessary in order that a w ao + a\ to the desired accuracy (i.e., up 
to error terms that are much smaller than aoP\/ pa\), and the former is certainly needed since 
E (N)/N < Anpao (i.e., the right side of (1.5) must be equal to Airpao plus a negative correction, 

which requires |oi| 3> ao\/ P a o)- 

For simplicity, we shall restrict our attention to interaction potentials given by the exponential 
function. Our proof can be adapted to a larger class of potentials, including the Yukawa 
potential. In our proof, however, we need the potential to be positive definite, with a Fourier 
transform satisfying nice decay properties as |k| — ► oo (e.g., polynomial decay), and our proof 
does not immediately extend beyond this class. Such restrictive condition is not supposed to 
have any physical relevance, of course, and (1.5) should hold for much more general repulsive 
potentials. We hope that the technical restrictions under which we proved Theorem 1 will be 
eliminated in future works [8]. 

The proof of Theorem 1 will proceed in two steps: we will get upper and lower bounds with 
the correct asymptotic form. The proof of the upper bound is based on a computation of 
the variational energy corresponding to the Bogoliubov trial wave function, following ideas of 
Girardeau and Arnowitt [5], see the next section. 

The strategy of the proof of a lower bound will follow closely the one of Lieb and Solovej in 
[6], where the ground state energy of bosonic jellium was investigated. We shall first localize 
the Hamiltonian in boxes of size £. Using the positivity of the Fourier transform of the ex- 
ponential interaction, we shall derive a preliminary estimate on the ground state energy and, 
correspondingly, on the degree of condensation in the small boxes. With this a priori bound 
on the number of particles n + outside the condensate, we shall be able to bound from below 
the full Hamiltonian by the Bogoliubov Hamiltonian minus an error term, depending on the a 
priori bound on n + . The key point is that it is possible to find a scaling regime for a and R 
such that the new error term is much smaller than the preliminary one, as Y — > 0. With this 
improved bound on the ground state energy we shall obtain new improved bounds on the size 
of fluctuations of n + that, in combination with the bounds for the ground state energy, will 
allow us to conclude the desired lower bound. 

II. THE UPPER BOUND 

Let us first derive an upper bound to the ground state energy, asymptotically agreeing with 
(1.5). In second quantized form, the Hamiltonian Hn can be rewritten as: 

H N = k 2 4c k + — - Y KP)4+ P c q-P c k c q , (2-1) 

k 1 1 k : q,p 
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where the sums run over vectors of the form lixvjL, v E 1? , and c k ,Ck are standard bosonic 
creation and annihilation operators, associated with the canonical basis of plane waves (for an 
introduction, see, e.g., [7]). Following [5], we choose the following variational state, inspired by 
Bogoliubov's approximate treatment of the weak coupling Bose gas: 

\SI b ,n) -cxp{i^^(k)(/3 1 a k -/3o4)}l^) ( 2 - 2 ) 

where: 

1) \CIn) — (A r !)~ 1 / 2 (co) Ar |0) is the ground state for N non-interacting particles; 

2) the operator a k is the pair annihilation operator a k = c k c_ k ; 

3) if we denote by iV = CqC the number operator associated to the constant wave function, 
Po is the partial isometry defined by 

SiT = coN ^ , = (/^ ) f = N~ V»4 , (2.3) 

having the properties 

f3 Q \n N ) = \si N _ 2 ) {n > 2) , pl\n N ) = ^Qn) = \n N+2 ) , 

[/3o,ATo]=2/3 , [(3o\N Q ]=-2fc 1 , 

[/?o,Ck] = [/?o,c k ]=0 (k^O); (2.4) 

4) ip is a continuous function from R 3 to R. 

Note that \SIb,n) is normalized, and that the particle number is equal to N. The variational 
principle implies the upper bound 

Eo(N) < (Ct BiN \H N \n BiN ) . (2.5) 

Following [5], after a lengthy but straightforward computation, we find that in the thermody- 
namic limit 



lim ^-(n B . N \H N \n B . N ) = \pv{0) / -^r[k 2 +p ^(k) + i/ 2 (k)l sinh 2 ?A(k) 

l^oo Jy 2 J (2ir)° L 2 J 

/ (03 [ P0iy(k) " ^ /l(k) l sinh ^( k ) cosh ^( k ) . ( 2 - 6 ) 



-p- 1 



where: 



p0 = p ~ j (03 sinh2 ^(q) 

7l ^ = / J0j* 1/ ^ k ~ ^ sinh cosh ^(q) 

/2(k) = / (0^ (k " q) sinh2 ^ (q) • (2J) 



Choosing ^>(k) = \ tanh 1 k^+^k) > we nn d that 



sinh^ ( k) = l k2+ ^-^ k4 + 2 ^ k2 , 
V 7 2 ^k 4 + 2/^(k)k 2 

Recall that f (k) is given in (1.4), and that a/ Rq = C^Y" 1 / 2 ^). A simple calculation shows that 

po = p(l + 0(y/^)) (2.9) 
for any d > 0. Moreover, we have the bounds 

|Ji(k)| < Cpa ^- . , |/ 2 (k)| < Cpa \[p~4 j (2.10) 

Ro [1 + (ki? ) 2 ] 2 V [l + (ki? ) 2 ] 2 
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for a suitable constant C. Substituting these bounds into (2.6) we find that 

Jim ^- T (n B , N \H N \n B , N ) = (2.11) 



= 1 



i- / 7^3 (k 2 + Mk) - v/k 4 + 2^(k)k 2 ) + o(v^) 

P 0,0 J R 3 (27T) 6 \ ) 



2np 2 

for < d < 1/4. A computation [7] shows that for d > the integral on the right side equals 

^ + -1^^3 + o( v ^). (2.12) 

Noting that a/(a + ai) = O(a /R ) 2 < F 1/2 for d < 1/4 this yields the desired result. 

Remark. The upper bound we have just derived yields the desired expression for any 
< d < 1/4. By suitably modifying the trial function ^>(k) above, one can actually show that 
the upper bound holds for any < d < 1/2 [3]. For d = 1/2, however, the ansatz (2.2) can not 
be expected to yield the Lee-Huang- Yang formula, even for the optimal choice of ip. 



III. THE LOWER BOUND 

We shall split the lower bound into several parts. The strategy is similar to the proof of the 
lower bound on the ground state energy of jellium by Lieb and Solovej in [6], and we shall refer 
to their paper for several essential ingredients. 

A. Sliding and localizing 

We start by rewriting (1.1) in the form 

N 

H' N = H N -4irNpa = -J2 A ^ + ^[ E ~ - (3-1) 



i=l u \<i<j<N 

N 



-pY] / dyv Ro (xi -y) + — // dx dy v Ro (x - y) 

,_, J A z J J Ax A 



with p — JV/|A|. We shall use the sliding method of [4] to reduce the problem to a small box. 

Let t, with < t < 1/2, be a parameter which we shall choose later to depend on p in 
such a way that t -> as pijg — »■ oo. Let \ € C(f(IR 3 ) be supported in [(-1 + t)/2, (1 - t)/2] 3 , 
< X < 1, with x(x) = 1 for x in the smaller box [(-1 + 2t)/2, (1 - 2t)/2] 3 , and x(x) = x(-x). 
Assume that all m-th order derivatives of \ are bounded by C m t~ m , where the constants C m 
depend only on m and are, in particular, independent of t. If M € N and I = M^ 1 L, let X£( x ) 
be a function on the torus A defined by X^( x ) = X)nez 3 x(^ _1 ( x + n ^))- For given x we also 
define 7 > by -f^ 1 = Jx(y) 2 dy, and note that 1 < 7 < (1 — 2t)~ 3 . We shall prove the 
following. 

Lemma III.l Let itR^ 1 be large enough. There exists a function of the form uj(t) = 
const t~ 1 Ro/£ such that if we set R^ 1 = R^ 1 +cj(t)/£ and 

«4(x, y) - X£(x)«fl(x - y)xt(y) ( 3 - 2 ) 
then the potential energy satisfies 

N p 2 r r 

Y v Ro (x t - Xj ) - pY / dyv Ro (x t -y) + // dx dy v Ro (x - y) > (3.3) 

l<i<j<JV 1=1 Ja J J Ax a 

-~w~ I 7J i wj?.( x i + z > x j + z ) - / d y w i( x j + z >y + z ) + 



me 



i£[l,...,M]3 Q m l<i<j<N j=l JA 

1 2 /•/■ a ,„,__. ,„a 

/Ax A 



+ ip 2 // dx dy (x + z, y + z) } - - const N^e^'^ , 
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where Q m is a cube of side length £ and centered at m£ (so that the collection {Q m }me[i,...,M] 3 
paves the torus A). 

The proof of Lemma III. 1 utilizes the following lemma, whose proof will be given after the 
proof of Lemma III. 1 . An analogous result for the Yukawa interaction potential was proved in 
[4, Lemma 2.1]. 

Lemma III.2 Let K : R 3 -> R be given by 

e -w|z| \ 

-—^) < M > 

with v > lo > 0. Let h satisfy (i) h is a C 6 function of compact support; (ii) h(0) = 1; (iii) all 
its m-th order derivatives, 1 < m < 6, are bounded by Ct x ~ m for some constants C > and 
t > 0. Assume further that h(z) = h(—z) so that K has a real Fourier transform. There exists 
a constant C\ (depending only on C but not on t, u or v) such that, i/min{l, ui}vt > C\, then 
K has a positive Fourier transform. 




Proof of Lemma III. 1 . We calculate 



U E / |^(x + z,y + z)= (3.5) 



Ro ^ Jo P 

U me[L..,M] 3 Jt "' 



~R^> Ja Xl<(K + Z ' ) VR ^ X ~ ^ + ^ = ~^o ht ^ X ~ ^ VR ^ ~ y ' ) ' 



/A 

-3 



where we have set he = j £ 3 xe * Xe- Note that he(x) vanishes if ||x|| > £, so we can naturally 
introduce a function h : R — > R of compact support and vanishing outside the cube of side 
2 centered at the origin, such that he(x) — J2nez 3 ^(^ _1 ( x + n ^))- Note that: (i) h(0) = 1; 
(ii) h has a quadratic maximum at z = 0; (iii) h is an even C°° function of compact support; 
(iv) all m-th order derivatives of h, m > 1, are bounded by C m t x ~ m , where the constants 
C m depend only on m and are, in particular, independent of t. The function h thus satisfies 
all the hypothesis of Lemma III. 2. Note also that the role of v and u> in Lemma III. 2 are 
here played by £Rq 1 and Lo(t) respectively. So, if w(t) > CiRq^ 1 ^ 1 , where C\ is the constant 
appearing in the statement of Lemma III. 2, we then conclude from it that the Fourier transform 
of the function K(x) = e~\ x ^ R ° - h(£- 1 x)e^ x ^ R (R/ Ro), is positive. Now, defining K A (x) = 
Snez 3 K(x + nL) and <p(x) = vr (x) — he(x)vn(x.)(R/Ro), we note that ip(x) = K A (x) + R(x), 
with |ii(x)| < const e~^ L ~^/ R . Because of positivity of the Fourier transform of K, 

E ¥>( x i ~ *j) ~ P^2 ^(x l -y) + -p 2 // ^(x-y)dxdy > 
i<i<j<JV j=i jA JJ « 3xl3 

> -^-K A (0) - const N* e -V-e)/R . (3.6) 

Since K A (0) < Ruj/£ + const cxp(-L/R ) this implies (3.4). ■ 

PROOF of Lemma III. 2. We write h(z) = l + q(z) + F(z), where q(z) is an even polynomial 
of degree 4 that vanishes at the origin, and F{z) < Ct~ 5 \z\ 6 . The Fourier transform of e — 
e -(i>+w)\*\ /(! + u / v } is given by 

8nu %TTV 4&TTl> 2 LJ 

+ p 2 ) 2 " ((!/ + UY + p 2 ) 2 - ((!/ + UJY + p 2 ) 3 ' 

Moreover, the Fourier transform of q(z)e~^ v+w ^ is 



(3.8) 
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whose absolute value, if vt > C\, can be bounded above by const • vt~ 1 [(i' + ui) 2 + p 2 ] -3 
(here we used that q is assumed to be even and that its m'th order coefficients arc bounded 
by Ct 1_m ). Finally we claim that the Fourier transform of F(z)e~(" +UJ ^ z \ is bounded by 
const • v~H-%u + lo) 2 +p 2 ]- 3 . To sec this, note that F (z.)er {v+u1 ^ is a C 6 function, and 
hence 

((1/ + uj) 2 + p 2 ) 3 J F(z)e-^+^e-^ z dz = J [((v + uj f - A) 3 F(z)e~^ +UJ ^ 

(3.9) 

It is not difficult to see that the latter integral is bounded by Ct 5 f 3 . After collecting all the 
terms, we arrive at the statement of the lemma. ■ 

Below, we shall choose the parameters t and I as functions of p, Rq and do- We shall choose 
them in such a way that t -C 1 and £ 3> Rq- Moreover, we will have conditions of the form 

ir^ 00 ' r^^ ' and w^ 1 (3 - 10) 

-Ko -K 1 paoV P% n a 

as pa^ — > 0, such that the error in the specific ground state energy corresponding to the term 
Nu>(t)R/(2l) in (3.4) is much smaller than Npaoy^ paf,, which is the precision to which we want 
to compute the ground state energy. 

Consider now the n-particle Hamiltonian 

I>gL + ^V, (3.1D 

3=1 

where we have introduced the Neumann Laplacian Aq^ z in the cube Q m , z = Q m + z and the 
potential 

W z (xi,...,x„) = ^2 «4( x » +z,x j +z) - / dy w% (xj + z, y + z) 

l<i<j<n j=i ^ A 

+ Jp 2 // dxdy^(x + z,y + z) . (3.12) 



2' 



Ax A 



Denoting by E^ z the ground state energy of the Hamiltonian if£, z in (3.11) considered as a 
bosonic Hamiltonian for n particles confined to the cube Q m ,z, and using Lemma III. 1 , we find 
that the ground state energy Eq of (1.1) can be bounded below by 

m£[l,.,M]3 ""*m ~ ~ 

Note that all the Hamiltonians z are unitarily equivalent to 

_^ A 0) + 7^r £ w ^( Xi , Xj )-pf2[ dy W *( Xj ,y) + lp> [[ dxdy^(x,y) 

(3.14) 

where A^ denotes the Neumann Laplacian for the j-th particle in the cube [— £/2,£/2] 3 . As a 
consequence, in the L — > oo limit, we have reduced the problem to studying the Hamiltonians 
R r l on £ 2 (H72,£,2] 3 "), given by 

n 

H? = -"£4 j) (3-15) 
+ E w R( x ^ x o)-py2 d y w R( x j,y) + oP 2 // rfxrfyw fl (x,y) 

^0 l<i<j<n j=l' / * 3 2 ^^M3x K 3 



7 



with w R (x,y) = x(x^- 1 )e"l x - y l' /fl x(y-0- If s " is thc ground state energy of ff™, from (3.13) 
we infer that 

^^« + ^"-^ <" 6 > 

In the remainder of this paper we shall study the Hamiltonians (3.15). For future reference, 
let us finally note that in second quantized form H™ can be rewritten as 

H'l = ^2 P 2fl p a P + ~^Jj4~ [\ Y ^pq,^ a P a q a M a ^ - P £3 *Op,o q apa q + ^p 2 i 6 w 00 ^ 
p pq,A*^ pq 

(3.17) 

where the sums over the momenta run over the values p = (j>i,P2,P3) such that ipi/n G 
Z + , a p , a p are bosonic creation/annihilation operators corresponding to the orthonormal basis 
</>p(x) = c-pfr 3 / 2 rij=i cos {Pj' K £~ 1 (xj + £/%)), and the coefficients w pq , M „ are defined as 

Wpq^u = J J rfxdyw fl (x,y)^ p (x)^ q (y)0 M (x)(/. I/ (y) . (3.18) 



B. Scalings 

Before proceeding with the proof of the lower bound, let us make a few remarks on the choice 
of the parameters ao,Ro, t, t. We recall t hat o ur purpose is to compute the ground state energy 
of (1.1) up to terms of the order Npao^/pa^, asymptotically as Y = pa% — » 0. In the following 
we shall choose a /Ro ~ y 1 / 2 -^ aQ /£ ^ yb+1/2 an j ^ ^ y T 5 where rf, b and t are positive 
scaling exponents. [Here by / ~ g we mean that C~ 1 g < f < Cg, for some universal constant 
C] We shall require d < 1/4. Note that the conditions b,d> and d < 1/4 imply in particular 
that ao <C i£o ^ an d a o /^ *S \/ P a o (two conditions that are of course necessary to be able 
to neglect finite size effects due to the boxes of size £) and that ao/ Ro 2> \J pa\ ^> (ao/Ro) 2 (a 
condition that is necessary for the Bogoliubov approximation to be valid, as explained in the 
Introduction). In order to be able to prove that the various error terms in our estimates are 
much smaller than N pao\/ pa^ we will be forced to require that b, d, r are small enough and that 
they satisfy a number of inequalities, some of which will now be discussed. Such inequalities 
will be satisfied by proper choices of b and r, as long as d is small enough. 

1. We require pRo » 1, that is < d < 1/6. Note that under this condition we also have 
ao/Ro > \fpal > (<V-Ro) 3 - 

2. We require URq 1 > 1, as in (3.10), so that b + d > r. 

3. We require aoRo 1 £~ 2 t~ 1 <C pa$y/ pa$, as in (3.10), so that 2b — d — t > 0. 

4. Noting that the contribution from the potential energy per particle in ff™ is expected 
(on the basis of Bogoliubov theory) to be of order nao-R^ 1 relative to the main term, we 
require both that tooi?^" 1 and ao(.R _1 — Rq 1 ) are much smaller than i/paf, in order to 
guarantee that the errors produced by the presence of 7 in front of the potential energy 
and by the replacement of i? with R are negligible. These conditions imply r > d and 
2b + d> t. 

Further requirements will be discussed below. 



C. A priori bounds on n and n+ 



As a first step in our argument, we shall derive preliminary bounds on the number of particles 
minimizing and on the average number of particles (n+) outside the condensate. [Here the 
operator n + is defined, in second quantized form, as n + = X) P ^o a p a p-] First of all, note that 
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m > 



> infi<„<Ar E™, so we can restrict our attention to the values of n such that E" < 0. As 
proved in the next lemma, such values of n cannot be too small, namely they all satisfy n > cp£ 3 
for a suitable constant c. We can thus assume, without loss of generality, that n > cpi 3 in the 
following. 

Lemma III. 3 If £/R and t^ 1 are large enough, then H™ > if n < p£ 3 /4. 
Proof. From the definition of Hf we see immediately that 

7a0 -[-pE / dyw R ( Xt ,y) + £ JJ d X dyw R (x,y)} > 2^<^- ( - 4np + p 2 £ 3 ) , 

(3.19) 

where we used that sup x J ior(x, y) dy < 8irR 3 and that, for i/R and t 1 large enough, 
Jj dxdywii(x.,y) > 4ttR 3 £ 3 . This proves the lemma. ■ 

A similar argument allows us to get a preliminary bound on the average number of particles 
outside the condensate (n+). 

Lemma III. 4 If t and (i?o — R)R^ 1 are small enough, then for any state such that (H™) < 0, 
the expectation of the number of excited particles satisfies (n+) < const nao£ 2 R^ 3 ■ 

Proof. Using the fact that the potential e~^/ R is positive definite, we obtain 

fl? > "5>?> -SE^.x.) • (3.20) 

The claim of the lemma follows by using 

(-ELi^ ) > {fi + )TT 2 /e and the fact that 

iojj(x,x) < 1. ■ 

Of course, in order for the bound in Lemma III. 4 to be useful, it must be oo£ 2 Rq 3 <C 1. In 
the following we shall impose this condition by requiring that, in terms of the scaling exponents 
introduced in the previous section, 2b + 3d < 1/2. We shall define vq = 1/2 — 2b — 3d, so that 
our preliminary a priori bound reads (n + )/n < Y v ° . 



D. Bound on the unimportant part of the Hamiltonian 

Motivated by Bogoliubov's computation of the ground state energy, we would like to be able 
to neglect in all terms but those containing precisely two a*, with p / 0. Let no = a ao 
and let us rewrite (3.17) in the form 



ja R 

p q 

ja R r 



H" = ^ p 2 a p a p + Y ^pq,oo (ala^aoao + 2a p a a a q + a^Opa^ + 



+ -^et |™oo,oo [(no - P^f - «o] + 2 ^ ™po,oo ("o - pP)a\>a Q + a o a p (n - p£ 3 ) 
p/o 



+ 



+ 2 Y w p o !q0 a p a q (no - P^ 3 ) + 2 ^ ^pq,f*o ^a P a q a M a + a a jL a q a p) + 

p,q/0 p,q,/x/0 
p,q,/x,^/0 

We would like to show that all terms but those in the first line are negligible. Let us then 
estimate all these contributions in terms of fi + and n — p£ 3 . We shall use, without proof, a 
number of lemmas from [6] . Note that from now on we shall always assume valid the conditions 
discussed in Sec. IIIB above. Note also that h = n — n + and |u) P q !/x „| < const R 3 /£ 3 . 
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1. The first term in the second line of (3.21) satisfies 

-^^woo.oo L( n o - P« ) -n J > const— (n -pf 3 ) -const— (3.22) 
and, for any £ > 

-^rWoo.oo [{no - pi ) -n \> -^^w 00i00 (l-e){n-pr) - const — ^ - const 

(3.23) 

2. By Lemma 5.5 of [6], the second term in the second line of (3.21) satisfies for any s > 

^S^^u> p0 ,oo (n - pe 3 )ala + ala p (h - p£ 3 ) (3.24) 

> —const — — const (no — pi 3 ) 2 — const e-^ 
^ w P o,oo [(no - p^ 3 )apa + a o a p (n - p£ 3 ) (3.25) 



and 



jaoR 



> — 51 ^pO,oo (n - /rf 3 )a p a + a a p (n - p£ 3 ) 
p#o 



„ a n n + a 

— const en+ — ^ const — - 

t 6 e t A 



3. By Lemma 5.3 of [6] the first term in the third line of (3.21) satisfies 

"^T" X! ^po,qoa P a q (wo - P^ 3 ) > -const ^[/rf 3 - n]+n+ - const °^ 3 + (3.26) 
p,q#o 

where [t]+ = max{t, 0}. 

4. By Lemma 5.6 of [6] and its proof, the second term in the third line of (3.21) satisfies for 
any e > 

■^p- X! ^pq.^o (a* a* a^ao + a a jx a q a P ) (3.27) 

p,q,/x/0 

a n „ 1 a n + 1 -fa R - t t 

> -COnst £ — n + - COnSt - — ^ W pq^^ a p a q a (i a " 

P,q,M,^/0 

5. The term in the fourth line of (3.21) satisfies 

"2^ ^ w pq ,^a p a q a Mai/ < const —^f- . (3.28) 

P,q,M:f#0 

This follows immediately from the fact that ior(x, y) < 1. 



Remark. According to Bogoliubov's theory we expect that in the ground state (n + ) <~ 
ny^paQ. From the upper bound in (3.28) we thus expect that the contribution to the ground 
state energy from the quartic term -j^- X) p q ^ ^pq,/ii> a P a q a /x a i> is at most ~ n 2 pciQ^. In 

order to show that Bogoliubov theory is asymptotically correct up to terms of order npao \J pa 6 ) 
we shall require such a bound to be much smaller than npa,Q\J pa^. For n <~ p£ 3 and in terms 
of the scaling exponents introduced above, this implies Y 1 / 2 ~ 3b ~ 3d ^ ^ ^sX is 36 + 3d < 1/2. 
In the following we shall assume this condition valid. It will be convenient to summarize here 
all the requirement we asked for so far on the scaling exponents introduced in Sec. IIIB: 

26 - d > t > d , and \>b + d>r . (3.29) 

6 

From now on we shall always assume that these relations are valid and that Y is small enough. 



10 



E. The quadratic Hamiltonian 

In this section we consider the main part of the Hamiltonian. This is the "quadratic" Hamil- 
tonian considered by Bogoliubov. It consists of the kinetic energy and all the terms with the 
coefficients w pq , o, Woo, Pq w p o,o q , and w 0p , q o with p,q^0, i.e., 

H B = ^2 P 2fl P a P + E *pq.°° ( a p a q a o»o + 2a^aJ,aoa q + <4<4a p a q ) . (3.30) 

p p,q^O 

In order to compute all the bounds we find it necessary to include the first term in the second 
line of (3.25) into the "quadratic" Hamiltonian. We therefore define 



tt r 2 t , 7 a o-R \ - „ 
Hq = p a pOp + —^r w p°°° 



+ 



p 

2m — 

p,q/o 



(n- pt i )aLa + a\a p {n - p£ 3 ) (3.31) 



^ w Pq ,oo (apa q aoao + 2a p a|,aoa q + alala p aqj . 



Note that Hb = Hq in the neutral case n — p£ 3 . Our goal is to give a lower bound on the 
ground state energy of the Hamiltonian Hq . 

For any keR denote Xe,k( x ) = e lkx x(x/£) and define the operators 

b\. = ^(0 p ,Xi?,k)apao and 6 k = ^ (Xl,k, M a o a p • ( 3 - 32 ) 
p#o P #o 

Note that they satisfy the commutation relations 

[bk, bl,} = n (xi,k, XLk>) - ^> k ' X^,k') a q a P - no(xe,k, <Po){<f>o,Xe,k>) ■ (3.33) 

p,q#o 

Using Lemma 6.2 of [6] we find that 



) 3 |k| 2 + (£t 3 ) 



(btbk) 



12 \"k 



(3.34) 



for a suitable constant C and for all states with particle number equal to n. Concerning the 
potential energy terms, note that we may write 



w fl (x,y) 



dk 



VR(k)x£,k(x)X£. k (y) . 



(3.35) 



where Vjj(k) = 87ri? 3 [l + (ki?) 2 ] 2 . The last two sums in the Hamiltonian (3.31) can therefore 
be written as 



^3 Vk(k) [(n - p* 3 )* 3 ' 2 (x(kO&l + X (W)fck) 



(2 



(3.36) 



+ 



k 6_ k + &£&L k + fe k fe_ k 



Rn 



E 

p,q/o 



t. 



Wpq,ooa p aq 



Thus, we have for states with particle number equal to n that 

ia R 



(Hq) > f (ft (k)) - ^ E ^ Pq ,oo (a p a q ) , (3.37) 

•/K3 l^TTJ «0 pq/0 



where 



M»0 = 



(1 - C'i) 2 



(fo k 6 k + fel k 6- k ) 



2n |k| 2 + (ft 3 )" 2 
+ 2W ^ (k) [ ( " " ^ 3)f3/2 (*M<& + + ** + fol k )) (3.38) 



+fo k fo k + &! k & -k + & k &t -k + fokfo-k 
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In order to give a lower bound on /iQ(k), we can use Bogoliubov's method, in form of Theorem 
6.3 of [6]. This theorem states that, for arbitrary constants A > B > and k £ C, the inequality 

A(blb k + fe f _ k 6- k ) + B(blbl u + fe k fe_ k ) + n(bl + fo_ k ) + k* (6k + fe f _ k ) (3.39) 



> -±(,4- v / ^ r ^)([ 6k , 6k ] + [ 6 _ k , 6 t_ k]) _ 



2k 



2 



holds. Note that in our case 



[&k, &L] - "° / dx x(x/^) 2 < n* 3 • (3.40) 



With the notation 



7 a i? - 
Bk = 2i^ (k) 



we thus obtain that on the subspace of n particles 



(l-C'tf |k| 4 
" 4k= 2n MHPP +Bk ' (3 ' 41) 

«k=2^2^«( k )(»-^)x(k0 



Moreover, since 



^2 ^pq,oo4 fl q = / f dyw R (x,y)[^2 ^W^fE ^p(y) a P < ( 3 - 43 ) 

!>-•<> p » 



p,q#0 p^O p/O 

< 



87Ti? 3 



we have that 

~^T~ ^pq,ooapa q < const . (3.44) 



p,q#o 



Using (3.37), (3.42) and (3.44), we find that, on the subspace with n particles, 

*> a -L ^ ( M - J* 1 *) + ^} - const "? ■ (3 - 45) 



Now, using Au > 2? k and the definitions of £> k , K k , we get 

|2 



/■ dk 2\ Kk \ 2 r jk_K^ 

As (2tt)3 A + £ k - J R s (2tt)3 S k 1 ' J 

7«0-R, ^2 /* dk ,>, n |M 1w 2 7«0-R, „3n2 - 

= W (n ~ p } As w * ( )lx( }l = W ( " " p } ™ 00 ' 00 ' 



As a result, on the subspace with n particles, 



H Q > -nl - "T^r (" ~ P^ 3 ) 2 ^oo,oo - const , (3.47) 



where 

1 r dk 

2^7 Wf 



' = ~ J J^f [ /(k) " ^f^Wf] , (3-48) 



f(k) = (1 - C>t?^ - „ + ^V R (k) 



n |k| 2 + (^ 3 )- 2 i?4 

Rn 
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Similarly, the Bogoliubov Hamiltonian in (3.30) on the subspace with n particles admits the 
lower bound 



Hb > —nl — const . 

t 6 



(3.49) 



Note that / > g > implies / — \J ' f 2 — g 2 < mm{g,g 2 /(f — g)}. Thus clearly / can be 
bounded as I < (2ir)~ 3 (2p)~ 1 / dkg(k) < const a R~ 3 , so that 



Hb > — const n-^ . 

R d 



Moreover, if n < Cp£ 3 , we find 



(3.50) 



|k| 2 <a p 



dk 

(2^F 5(k) 



dk 



;(k) s 



< const i pa J pa 3 , + pa\ \ dk 



|kp>a„ P (2^) 3 /(k) - ff(k) 

1 

,-„„ l(kR) 2 + I] 4 



(1 + (Ht 3 ) 



< const pa i J pal + ^ 



R 



1 



/pa^R V^ 3 



If the scaling exponents satisfy 



2b + d- 6r > , 



(3.51) 



(3.52) 



then the last expression in (3.51) can be bounded from above by const pa ^- Hence, if n < Cp£ 2 
and 26 + d - 6r > 0, 



Hb > — const npa . 



(3.53) 



F. Improved bounds on n 

Using the bounds derived in the previous sections, we shall now get an improved bound on 
n, which implies that for states with (H™) < 0, n cannot deviate too much from p£ 3 . In order 
to bound (3.21) from below, we use (3.22), (3.24), (3.26), (3.27) and (3.50) [note that we shall 
use (3.27) with e replaced by £ -1 ]. The result is that, for some positive constants c and C, we 
have 

H?>(c-Ce)^(n -p£ 3 ) 2 + (l-Ce)^ ]T » Pq , F «f<fl,»„ (3.54) 
f na n + na , a a 3 „ a n+ a ri + a 

+ T— + e F + ^ n+ + ^ + e -ffr + n ^ 

for some e > 0. Choosing e = min{l, c}/(2C), using n + < n, E P . q .^.i/^o ! "pq^^ fl p a q fl /i a " ^ 
and recalling that, by Lemma III. 3, n < p£ 3 /4 implies H™ > 0, we have, for some new constants 
c' and C and for any state with (if™) < 0, 

> m > c'^{(h - P e?) - C'^ + } (3.55) 

and, therefore, 

(n-p£ 3 ) 2 f<n+) 1 1 ,„ r „x 

< const <^ ^± + 3.56 

[_ n pR A J 
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Here, we used ((no — p£ 3 ) 2 ) > {n — pi 3 ) 2 — 2n(h + ). Now, let us recall from Section IIIC that, 
in terms of the scaling exponents b, d, we have (n + )/n < yi _2l)_M , and (pR 3 )^ 1 ~ Y?~ 3d , so 
that 

( "~f 3)2 < const Y va , (3.57) 
n^ 

where = 1/2 — 26 — 3d as before. Eq. (3.57) can be rewritten as 

\n - p£ 3 \ < const pi 3 Y Ua/2 . (3.58) 

In order to get the bounds above we sacrificed all the kinetic energy in (3.21). Of course 
this is not necessary: we can decide to sacrifice only half of it and we would still get the same 
bounds, only with different constants. If we proceed in this way we see that for any n-particle 
state such that (Hf) < 0, 

^2 |p| 2 (« P a P> < const npa Y v ° . (3.59) 



G. Localization of n+ 

The idea now is to use the improved bound on n together with the bounds in previous sections 
in order to find an improved bound on the energy of the ground state. In order to do this it 
is clear from the bounds in Sec. HID that we need to estimate (n+). Since we have bounded 
only n+ so far, we would like to argue that (n+) ~ (h + ) 2 . In this section we shall discuss how 
to do this. We shall utilize the following theorem, which is Theorem A.l of [6]. [The fc th supra- 
(resp. infra-) diagonal of a matrix A is the submatrix consisting of all elements a^i+fc (resp. 

Theorem III.l Suppose that A is an N x N Hermitean matrix and let A k , with k = 
0, 1, . . . ,N — 1, denote the matrix consisting of the k th supra- and infra-diagonal of A. Let 
ip G C N be a normalized vector and set d^ — (ip,A k 4>) and A = (tp,Aip) — X^fcLo 1 ^ k nee d 
not be an eigenvector of A). Choose some positive integer M < N. Then, with M fixed, there 
is some n G [0, N — M] and some normalized vector (f> G C N with the property that <f>j = 
unless n + l<j<n + M (i.e., 4> has length M) and such that 

„ M-l N-l 

(0,^)<A + — £ fc 2 |4|+C^ |4| (3.60) 

fe=l k=M 

where C > is a universal constant. (Note that the first sum starts with k = I). 

From this theorem we can get a localization bound on n + in the following way. Consider a 
normalized n-particle wavefunction 'J, which we may write as ^ = X)m=o Cm * m ' where for all 
m = 0, 1, 2, . . . , n, is a normalized cigenfunction of n + with eigenvalue m. We now consider 
the (n + 1) x (n + 1) Hermitean matrix A with matrix elements A m m' = ('I'm, H^ m i). 

We shall use Theorem III.l for this matrix and the vector ip = (cq, . . . , c n ). We shall choose 
M in Theorem III.l to be of the order of the upper bound on (n + ) derived in Lemma III. 4, 
e.g., M is the integer part of nY va . Note that, if n ~ p£ 3 , we have M ^> 1. With the notation 
in Theorem III.l we have A = (ip, Aip) = Hfity). Note also that because of the structure of 
H™ we have, again with the notation from Theorem III.l, that dk — if k > 3. We conclude 
from it that there exists a normalized wavefunction VP with the property that the corresponding 
n + values belong to an interval of length M ~ nY VQ and such that 

const —L-fldil + ldal). (3.61) 

We shall now bound d\ and ^2- We have d\ = ($f, H"(l)^>), where H™(1) is the part of the 
Hamiltonian H™ containing all the terms with the coefficients w pqi)11/ for which precisely one 
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or three indices are 0. These are the terms bounded in (3.24) and (3.27). These estimates are 
stated as one-sided bounds. It is however clear that they could have been stated as two sided 
bounds. Using in addition the bound (3.28) and h\ < nh + then, for any e > and some 
positive constant C, we get 

Mil < C(% e-^(no pi 3 ) 2 +s^+ e^±\ *) . (3.62) 

If $ satisfies (\I>, Hfty) < 0, we can use Lemma III. 4 and (3.58) to conclude that 

Optimizing over e > yields the bound 

K| < Cnpa n Y"°-'i {b+d) = Cnpa Y^-i b -'i d . (3.64) 

For d 2 we obtain 

\<k\ ^pq.oo (a^aoao + ao a o a P«q) 

p,q/o 

= [J2 P 2fl P a P + ^ u>Pa,ooay o a a^ tt) - (*, H B ^) (3.65) 



p,q/o 



where 

jaoR 



Hb = ^p 2 aj,a p + ^ r ^2 ^pq,oo (-<4,a q a a + 2aJ,aJ,a a q - alala p a^j (3.66) 

p p,q#0 

is an operator unitarily equivalent to Hb- (It is obtained from it by replacing a p ,a p by 

— ia p ,ia p , respectively.) Of course Hb satisfies the same lower bound (3.53) as Hb- It is 
not difficult to see that 

< -5^" V ^pq,oo(ai,ao a oaq) < — 73— 5l"(«+) (3.67) 
^° P ,q/o 1 H " 

(compare with Lemma 5.4 of [6]). If <J> satisfies ('J, H™^) < then, using (3.65), (3.67), Lemma 
III.4, (3.59) and (3.53), we get: 

\d 2 \ < const npa a jY" + ^} < const npa Y Va - (3.68) 

Putting together these bounds we find that if (*, £f™\E') < then there exists a normalized 
wavefunction "J" with the property that the corresponding n+ values belong to an interval of 
length M ~ nY Va and such that 

(*, H^) > (tf , Hg^j - const -^-F^-K b + d ) 

> i/; 1 *) - Cnpa r M0 , (3.69) 

where p = — ^0 + 1 + — |^ = ^ + ir^ + Since /x > 1/2, the error term in the 



last line is much smaller than npan\/ p%. Without loss of generality, we may assume that 
< 0, in which case Lemma III. 4 implies that (^>,h + $>^ < const nY v ° . We also 

know that the possible n + values of ^ range in an interval of length M ~ nY Va . This implies 
that if ( , i , ,i^™ , J , ) < —CnpaoY 1 / 2 then the allowed values of n + for 4> are less than CnY Va , 
for a suitable constant C. In particular, (n+) < Cn 2 Y 2v ° in the state \E r . Hence, as far as the 
derivation of a lower bound on the ground state energy is concerned, it is not a restriction to 
assume that {n\) ~ (n+) 2 ■ This fact will be used in the next section to derive improved lower 
bounds on the ground state energy. 
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H. Improved bound on the ground state energy 



Let \P be the ground state of H". In this section we shall get an improved lower bound on 
( l i',H^) under the assumption that (vE", -ff"\I/) is small enough such that (3.69) implies that 
(\I/, H™$>) < 0. Note that if this assumption is violated then the desired bound on the ground 
state energy would automatically be true. Hence, as discussed at the end of previous section 
we know that 



(*, fl?*) > (H?) - Cnpa 



(3.70) 



where the average (•) is over an n-particle state with allowed values of n + smaller than CnY l/ " . 
Using (3.21), (3.23), (3.25), (3.26), (3.27) [this time precisely in the form stated, without 
replacement of e by and (3.28), we find that 



(fl?) > (H Q ) + ^^oo,oo (1 - e)(n - pi 3 ) 2 - C 



2Rt 



«o («+) a n a n 



+^\n- P e 3 \(n + ) + 



£ 3 



aon 



1 a Q (n + ) 1 a (n+) 



i? 3 



+ e" R 3 



(3.71) 



Now, if < e < 1, using n + < CnY v ° , (3.47), and and (3.58), we get from the last inequality 
that 



EY Va + Y^ Va + _y 2i/ o — 36— 3d 

£ 



{Hi) > -nl - Cnpa 
Optimizing over e yields 

(#") > -nl - Cnpa Y ai 
6d. If 



where ot\ = | 



^_3 6 _3 d= 3_9 & . 



3 9, „ , 1 

4-2 6 - 6d> 2' 



the error term npaoY ai is much smaller than npaoY 1 ' 2 and, therefore, 

E% > -nl - npa Q o(Y^ 2 ) . 



(3.72) 



(3.73) 



(3.74) 



(3.75) 



We are left with estimating the constant I defined by (3.48). It is not difficult to see that, 
under the assumptions made so far on the scaling exponents, 



I = 47r j oa I — 
»o 



(3.76) 



The conditions (3.29), (3.52) and (3.74) on the scaling exponents that we required for the proof 
to work can be summarized into the following conditions: 



2b + d > 6r 



t > d , 



- > 36 + Ad . 
6 



(3.77) 



It is easy to check that if d < 1/69 then all these requirements on the scaling exponents can be 
satisfied. ■ 
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